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Abstract 

We give a detailed mathematical analysis of the radiative transport limit for the average 
phase space density of solutions of the Schrodinger equation with time dependent random 
potential. Our derivation is based on the construction of an approximate martingale for the 
random Wigner distribution. 

1 Introduction 

The Schrodinger equation with random potential arises in many applications, especially in wave 
propagation in random media, in the paraxial or parabolic approximation. In this case the 
time harmonic wave field has the form u = e*"^~"*(/)(z, a;) where k, = w/c is the free space 
wave number, z is the coordinate in the direction of propagation, x are the coordinates in the 
transverse directions and (j) satisfies the Schrodinger equation 

2iK^ + Aj;^ + = 0. (1) 

Here ij.{z,x) — n'^{z,x) — 1 denotes the fluctuations of the index of refraction. The original 
scattering problem for (A + K^n^)u — becomes an initial value problem for ip in z, in this 
approximation, so (j) must be given at z — 0. The validity of the parabolic approximation in 
random media under different scaling limits is considered in jl], ^ and in general in Q]. 

The purpose of this paper is to prove a theorem that establishes the validity of the transport 
approximation for the average Wigner distribution of </>, in a suitable scaling limit and for a class 
of stochastic models for the index of refraction fluctuations n that are Markovian in z. For this 
class of stochastic models it is possible to use martingale methods to simplify the analysis. 

Since the coordinate z in the direction of propagation plays the role of time in the parabolic 
approximation we will denote it by t in the rest of the paper. The problem then is to analyze the 
Schrodinger equation with time dependent potential and to show that the associated average 
Wigner distribution converges to the solution of a radiative transport equation. 

The propagation of wave energy in a scattering medium is described phenomenologically by 
radiative transport theory Q as follows. Multiple scattering creates waves with all wave vectors 
fc S R'^ at every position x £ S.'^. Let us denote by W{t,x,k) the energy density of a wave 
having wave vector k at position x at time t. The energy balance equation has the form 

^^^t'^''^^ +k-V^Wit,x,k) = I dp a{x,k,p)W{t,x,p) -j:{x,k)W{t,x,k). (2) 
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Here a{x, k,p) is the probability to scatter from direction p into direction k at position x, and 
k) is the total probability to scatter from direction k into some other direction. When 
energy is conserved and scattering is symmetric, 



'S{x,k)— / dp<7{x,k,p), (t{x, k,p) = (t{x,p, k) 



equation (Q) may be rewritten as 
dWit,x,k) 



dt 



+ k ■V^W{t,x,k) ^ dp cr{x,k,p)[W{t,x,p) -W{t,x,k)]. (3) 



While various formal derivations of the radiative transport equation (g), starting from the 
wave equation in a random medium, have been known since the mid-1960's (see [ p2| for an ex- 
tensive bibliography) the mathematical methodology for doing this is not very well developed. A 
rigorous derivation of a spatially homogeneous transport equation starting from the Schrodinger 



equation is given by H. Spohn |23 who derived (||) for sufficiently short times t € [0,fo] and 
with a time independent Gaussian potential. This result was extended to higher-order correla- 
tion functions by T. Ho, L. Landau and A. Wilkins with the same restrictions. Recently L. 
Erdos and H-T. Yau jl^ removed the small time restriction and considered more general initial 
data. The idea in these proofs is to consider the Neumann series expansion for the solution of 
the Schrodinger equation and to infer appropriate estimates from it that allow passage to the 
limit. 

In this paper we deal with time dependent random potentials for which it is possible to 
analyze the transport approximations in a relatively simple manner, without infinite Neumann 
expansions. We model the random potential by a Markov process in time so that we can use 
martingale methods and suitable test function expansions. A limit theorem for one dimensional 



waves where such methods are used is given in 1 20 and more general ones in [|6| . An analysis of 



the random Schrodinger equation with rapidly decorrelating in time potential is given in |21 
Limit theorems for linear random operator equations that decorrelate rapidly in time are given 
in JT9[ . The random Schrodinger equation with delta function potential is analyzed in Q using 
martingale methods and its equilibrium solutions are constructed in . 

Formal derivations of radiative transport equations for various types of waves in random 
media are given in ^ . Appendix |^ contains such a derivation for the time-dependent 

case. The results presented here extend to linear hyperbolic systems with random coefhcients 
such as those considered in [ p2| . The transport equation for the Schrodinger equation with a 
time-dependent potential was used recently in 0| to explain pulse stabilization in time-reversal 
in a paraxial approximation to the wave equation. The same phenomenon for general wave 
equations in time-independent media was related to transport theory in Full justification of 
the results in Q requires analysis of the higher-order correlation functions and specific scalings, 
as discussed in p^ . 

Acknov^rledgment. We thank K. Solna for numerous discussions. G. Bal was supported 
by NSF grant DMS-0072008, L. Ryzhik by NSF grant DMS-9971742 and G. Papanicolaou by 
grants AFOSR F49620-01-1-0465 and NSF-DMS-9971972. 

2 The main result 

2.1 The Wigner distribution and the main result 

We consider the initial value problem for the Schrodinger equation (Q) in dimensionless form 

+ 4a0, -V^vfi, -) 0, = 0. (4) 



dt 2 ' \ee 
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The initial data 4'^{x) = (/)e(0,x) are assumed to be uniformly bounded in L^(R'^): 

< C. (5) 

This implies that 

\\Mt)\\L^ = Ue\\L-<C, t>0 (6) 

since the i^-norm of the solution is preserved by the Schrodinger evolution (^. Conservation 
of the L^-norm implies that for every realization of the random potential Ve{t,x) = V (tj e,x j e) 
there is a sequence £fe ^ so that the energy density 

E,^{t,x)^\<l>,,(t,x)\^ 

has a weak limit E[t, x) in i^([0, T] x R'^) as fc oo. 

It is well known that the limit E(t, x) does not satisfy a closed equation. A convenient way 
to study this limit is to consider energy propagation in phase space that includes all positions x 
and wave vectors k using the Wigner distribution We(t, x, k), defined by 

w,{t,x,k)^[ J^e^'^-yMt,x-'-^mt,^ + '-^)- (7) 

Here * denotes complex conjugation. The basic properties of the Wigner distribution can be 
found in p3|, p^. In particular, 

dk We{t,x,k) = \(t>,{t,x)\^ = E^{t,x) 

but may not be interpreted as energy density in phase space since it is not necessarily 
positive. However, the limit of along a sub-sequence ^ exists in iS'(M'^ x ffi.''), the space 
of S chwartz distributions, and is positive |13, 16|. We recall the proof of convergence in Section 



2.3. We say that the family ^^(x) is pure if the family of its Wigner distributions We{x,k) 
converges weakly as £ ^ to a distribution Wo(x, k) € S' without restriction to a subsequence. 
We will assume that the initial data form a pure family throughout the paper. 

Our main result concerns the convergence of the expectation of the Wigner distribution 
We(t,x,k) defined by (0) to the solution of the radiative transport equation (||). 

Theorem 2.1 Let the random potential V{t,x) satisfy the assumptions in Section \2.^ below. 
Let Wo{x,k) be the limit Wigner measure of the family (/)^(x), and let W{t,x,k) be the weak 
probabilistic solution of the transport equation 

dW 

_ + fc • V^W = CW (8) 

with initial data Wo{x, k) and where the operator C is defined by 

C\ = j j^R{^^,P~k)i\ip)~m)- (9) 



Here R{lo,p) is the Fourier transform of the correlation function of V , defined by ^Ijj). Then 
the expectation E {W^{t, x, k)} of the Wigner distribution W^{t, x, k) of the family (j3g{t,x) of 
solutions of (0j converges weak-* in L°°([0, T];S'{W^ x R'')) to W{t, x, k). 

By weak probabilistic solution we mean that W{t, x, k) satisfies 

dX 



{W,X)it)-{Wo,X\t=o)= {w, — + k-V,X + CXys)ds (10) 

for all test functions A € C^([0, T];S). We actually show that E {We(t, x, k)} converges to W in 
a smaller space £°°([0, r];^'(R'^ x R'')). The space A', which is defined below in Section U, is 
more convenient for the analysis of the Wigner distribution than S' . 
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2.2 The random potential 

The random potential V{t, x) is assumed to be stationary in space and time and to have mean 
zero. It is constructed in Fourier space as follows. Let V be the set of measures of bounded total 
variation with support inside a ball Bj^ = {\p\ 5: L} 

V = |f : J ^ \dV\ < C\ supp V C Bl, V{p) = V*{-p)^ (11) 

and let V(t,p) be a mean-zero Markov process on V with generator Q. The time-dependent 
random potential V{t, x) is given by 

and is real and uniformly bounded: 

\V{t,x)\<C. 

We assume that the process V{t, x) is stationary in t and a: with correlation function R{t, x) 

E {V{t, x)V{t + s,x + z)} ^ R{s, z) for aU x, z e W^, and t,s£M.. 

In terms of the process V{t,p) this means that given any two bounded continuous functions (j}{p) 
and '0(p) we have 



E- 



[{V{s),i){V{t + s),i;)] = {2^Y j dpR{t,p)cj,{p)i;{~p). (12) 

Here (•, •) is the usual duality product on M'' x M'^, and the power spectrum R is the Fourier 
transform of x) in x: 

We assume that R{t,p) £ S{M. x R'^) for simplicity and define R{uj,p) as 



R{uj,p) = J dte-''^'Rit,p), (13) 

which is the space-time Fourier transform of R. 

We assume that the generator Q is a bounded operator on L°°(V) with a unique invariant 
measure tt(V) 

Q*7r = 0. 

and that there exists a > such that if {g, tt) — then 

lle'^Slk^ < qiglk^e^"*. (14) 

The simplest example of a generator with gap in the spectrum and invariant measure tt is a 
jump process on V where 

QaiV) = / g{Vi)d7r{Vi) - g{V), f dniV) = 1. 

Given (|lj), the Fredholm alternative holds for the Poisson equation 

Qf = 9, 

provided that g satisfies {TT,g) = 0. It has a unique solution / with (tt, /) = and ||/||l^ £ 
C||(/||ioo. The solution / is given explicitly by 



fiV) - - / dre'-^giV), 
Jo 

and the integral converges absolutely because of (|l4|). 
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2.3 General convergence of the Wigner distribution 

Existence of the limit of the Wigner family W^{t,x,k) defined by is shown as follows. We 
introduce the space A, as in fl^, of functions A(a;, k) of x and k such that 



1A|| 



A 



where 



dy sup |A(a;,?;)| < oo, 



dke-^^-y\{x,k) 



(15) 



is the Fourier transform of A in k. Convergence in the space A is easier to establish than in S 
because its definition does not involve derivatives. Moreover, as the following lemma shows, the 
distributions are uniformly bounded in A' , the dual space to A. 

Lemma 2.2 The family a;,fc) is uniformly hounded in A', that is, there exists a 

constant C > independent of t so that: 



m{t)\\A' <c 



for all e > and t > 0. 
Proof. Let A(a;, k) e A. Then, 

{WS)A)= J dxdkW,it,x,k)X{x,k)^ I ^gM^e*Vai,^-y) 

^^^(t)eit,X+YWeit,X- y)A(x,?/). 

Therefore, using the Cauchy- Schwartz inequality in x 



(16) 



l(t,x+'-^)Xix,k) 



\{WS),X)\< 



dydx 
(2^ 



X{,x,y) (j)eit,x- 



2 ' 



sy 

(j)e{t,X+ y) 



< 



1A|U<C||A|U, 



where we use the conservation of the L^-norm (||) in the last step. This gives (p^). 

Lemma 2.2 implies that at every time i > we can choose a sequence Sj ^ so that W^^ 
converges weakly in A' C 5' to a limit distribution W{t). One can show that the limit 
measure W(t) is non-negative and may thus be interpreted as the limit energy density in phase 
space. Moreover, if there are no oscillations in the initial data on scales smaller than e then the 
limit captures correctly the behavior of the energy £'e(i, a;). More precisely, if 

e\\Vc^"jL2<C (17) 



then for any test function 9{x) G 5(]R'') we have 

dxdk0{x)Wit,x,k) = lim J dx9{x)\<j)e{t, x)\'^ . 
Condition (bjh is sufficient but not necessary for this. 



3 Convergence of the expectation 



We prove Theorem |2.l| in this section. The proof is based on the method of [gO| and proceeds 
as follows. The distribution W^. defined by (0) satisfies 



k ■ V^Wr 



dW, , 1 

W,iO,x,k) = W°ix,k). 



dV{t/e,p) 



Weit,x,k-^)~We{t,x,k + ^) (18) 
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Here is the Wigner distribution of the family (x) , the initial data for . The Cauchy 
problem generates a measure on the space C{[0,T]; A') of continuous functions in time 
with values in A' . It is supported on paths inside a ball X — {W G A' : < C} with the 

constant C as in (p^. The set X is the state space for the random process We(t). The joint 
process {V{t/e), We{t)) takes values in the space VxX. We will denote by the corresponding 
measure on the space VxX generated by ( p^ ) and the process V{t/e). Let us fix a deterministic 
function A € Ci([0, T]; 5). We will show that the functional Gx : C([0, T];X)^ C[0, T] defined 

by 

Gx[W]it)^{W,\)it)~ I (W, — + k-V,\ + CX)is)ds 







is an approximate P^-martingale. More precisely, we will show that 

{Gx[W]{t)\Ts} - Gx[W]is)\ < Cx,tV^ (19) 

uniformly for all 1^ £ C([0, T]; X) and < s < < < T. 

The next step is to show that the measures Pg form a tight family, so that there exists a 
subsequence ej ^ so that Pg converges weakly to a measure P supported on C{[0,T];X). 
Weak convergence of P^ and the strong convergence (|l^) together imply that Ga[W^]{*) is a 
P-martingale so that 

P^ {Gx[W]it)\T,} - Gx[W]{s) = 0. (20) 

Taking s = in this we obtain the transport equation (||) for W = {W{t)}, in its weak 
formulation (p^). 

The limit measure P may depend on the choice of the subsequence Sj —^ but the expectation 
W being the unique solution of (||) does not depend on it. Therefore the whole family SjWe} 
converges to as e ^ in L°°{[0,T];S'). Furthermore, the a priori bound ( |l6|) implies that 
We{t,x,k) converges weak-* in L°°{[0,T];A') for every realization of the random potential. 
Therefore the result above implies that actually E {We} converges to W in L°°{\0, T];A'). 

The proof is organized as follows. The approximate martingale property (|l9|) is proved in 



Sections 3.1 and 3.2. The weak compactness of the family P^ is proved in Section 



3.1 Construction of the test functions 

In order to obtain the approximate martingale property (^) one has to consider conditional 
expectation of functions F{W, V). The only functions we will need to consider are those of the 
form F{W,V) = {W,X{V)) with A € L°°{V;C\[0,T];S{R'^ x M''))). Given a function F{W,V) 
let us define the conditional expectation 

e!'^., \ F(W, V) !• (t) = E^- \ F(W(t),V(t))\ W(t) = Vit) = F S- , t > t. 



w,v 



[f{W, V)] (r) - E^^ [f{W{t),V{t))\ W{t) = V{t) - v] , 



The weak form of the infinitesimal generator of the Markov process generated by P^ is given by 



and hence 



Gl - {W, \iV))it) - ^ (w, + ^ + fc . V, + -^^V, A^ is)ds 



is a Pg-martingale. The operator /C is defined by 

1 f dV{p) 



IC[V, z]^{x, z,k,V)^- I ^^e'^-^ i^ix, z,k^?-)- i,{x, z,k + ?-) . (22) 
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The generator ( |21| ) comes from equation ( |18[ ) written in the form 

dtW, + k ■ V,W, = -^IC[Vit/e),x/e]W,. 

Given a test function X{t, x, k) £ C^{[0, T];S) we construct a function 

Xs{t,x,k,V) = X{t,x,k) + ^/eXl{t,x,k,V) + eX^it, x, k,V) 



(23) 



(24) 



with Af 2 (t) bounded in L°° (V; A{R'' x M'')) imiformly in t £ [0, T] . It is sufficient for us to prove 
the simpler bound for the correctors in A instead of S because of the a priori bound ( p^ for 
in A' ■ The functions 2 wiU be chosen so that 

for ah f G [0, T]. The approximate martingale property ( [l9| ) follows from this. The approximate 
test function Xe{t,x,k) in (|2^ ) is constructed in a manner similar to the formal asymptotic 
expansion (^3|) considered in Appendix 

The functions X\ and A| are as follows. Let Ai(i, z, /c, V^) be the mean-zero solution of the 
Poisson equation 

■ V^Ai +QAi = -O. (25) 

It is given explicitly by 



Xi{t,x,z,k,V) = - dre'''^ 



(2nY ' 



X{t, X, /c - |) - X{t, fc + |) 



Then we let X\{t,x,k,V) — Xi{t,x,x/e,k,V). The second order corrector is A|(t, x, /c, F) 
X2{t, x, x/e, k, V) where A2(i, a;, z, k, V) is the mean-zero solution of 



V,A2+QA2 = £A-Oi, 
which exists because E {/CAi} — CX, and is given by 

X2{t,x,z,k,V) = - dre"^ CX{t, x,k) ~ [ICXi]{t, x, z + rk,k,V) 
Jo L 

Using (pa) and (Un) we have 



(26) 



d_ 

dt 
d_ 
dt 



h=0 

k-V^] X + CX 



= {W,{^^+k-W^ + -^nv, J + ig ) (A + V^Xl + eXl 



d_ 

dt 



k-VA {V^X\ + eXl) + ^elC[V, -]Xl 



+ k-\iAx + cx) + ^/?.{W, c^) 



with 



d 



dt 



d 



dt 



The terms k ■ Va;Af 2 above are understood as differentiation with respect to the slow variable x 
only, and not with respect to x/e. It follows that G\ is given by 



Glit) = {W{t),X,) 



ds(W,\ ^^+k-W^ + C]X){s) 



ds{WXe){s) 



and is a martingale with respect to the measure defined on D{[0,T]; X x V), the space of 
right-continuous paths with left-side limits j|] . The estimate (|9|) follows from the following two 
lemmas. 
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Lemma 3.1 Let A € C^{[0,T];S{'M.'^ x R'^)). Then there exists a constant Cx > independent 
of time t g [0,T] so that the correctors Af (i) and X2{t) satisfy the uniform bounds 



and 



dt 



\\m)\\L^(V;A) 

k-V^\l{t) 



dt 



(27) 
(28) 



Lemma 3.2 There exists a constant Cx such that 

\\lC[V,x/e]\\A^A<C 

for any V G V and all e E (0, 1] . 

Indeed, iM) implies that \{W,X) - {W,X^)\ < C^/e for all W e X and V e V, while (| 
and Lemma 3.2 imply that for all t e [0, T] 



KemA<C 

for aU y e V so that (Oh follows. 

Proof of Lemma 3.2. Let ge{x, k) = IC[V,x/e]ri{x, k) with ri{x, k) E A. Then 



(29) 



fj{x,y) 



g-ip-j//2 _ gip-J//2 



and thus 



\9e{x,v)\ < C\fiix,y)\ 



and the conclusion of Lemma 3.2 follows. 



3.2 Bounds on the correctors 

We now prove Lemma We will omit the time dependence of the test function A to simplify 
the notation. 

Proof of Lemma |3.lj. We first prove (p^) . The Fourier transform of Af in k is given by 



Xl{x,y,V) 



1 



dre^'^ I |^A(x,y-rp)e"/-^ 



^-ip-{y-rp)/2 _ ^ip-{y-rp)/2 



Therefore using (14) we obtain 



||Af(x,y,T>)|U^^^ <-||A|Uo 



(30) 



(31) 



uniformly for all y G V. It is therefore sufficient to consider \y\ > 2. Let S{y) = {\y\ — 1)/4L 
with L as in the definition dill) of the set V. We write (pQ) as 



Ki^' ^) — Jr<S{y) + Jr>S{v) 



with 



J, 



r<S{y) 



siv) 



dre 



^-ip-(y-rp)/2 _ ^ip-(y-rp)/2 
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and 



Jr>Siv) = - r dre^"^ f ^^K^,y - rp)e"''-'' 



s{y) 



(27r)" 



g-ip-(y-rp)/2 _ ^ip-(y-rp)/2 



We estimate each of these two terms separately. 

To bound J^^g^y-j we note that since A is of the Schwartz class we have for \p\ < L and 
r < Siy) 



\{x, y - rp) 


< sup 


A(a;, z) 


< sup 


A(a;, z) 


< 




\z-y\<rL 




NI>lal/2 







Then we obtain 



dy sup 

\v\>2 



Jr<S{y)ix,y,V) 



|5rf" 



-4 / dre-'^- < Cx,c,. (32) 



\y\>2 \y\ Jo 



Next we note that 



dysup|J,>5(,)(x,y,y)| < |1A|U» / dy / dre""'^ < C„||A|U=c^. (33) 

|y|>2 x,V ' "'|y|>2 Js{y) 



Therefore (31), (B2) and (B3) imply that 



(34) 



for all V eV. 

We show next that Al is uniformly bounded. This is done in several steps that we formulate 
as separate lemmas. Define 

I{x,y)^ Tdre-^'-sup / T dse-"^ sup f \^^^\~X{x,y - rp - (r + s)q)\. 

Jo V J v^"""/ "'O Vi J y^''^! 



Lemma 3.3 We have the estimate 

\\\(x,y,V)\<C^ /(x,y) + |£A(x,j/)| 



(35) 



Lemma 3.4 For the limit operator C we have the hound 

\\C\\\a < I dS,dy\Cr\{U)\<C I d^dy\\{U)\- 



Here f denotes the Fourier transform of f both in x and k. 
We split I{x, y) as 

jS<S(y) , jS>S(y) 



l(X,y) - lr<S{y) +U>S(y) - ^r<S{y) + ^r<S{y) +V>S(y)- 



(36) 



(37) 



Lemma 3.5 We have the following bounds: 



dy sup 



Ir>Siy)ix,y)+K<IU\<CM 



and 



dy sup /;<|g; < C„ sup I |2/nA(x, y)\} 



(38) 
(39) 
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Lemmas 3.3, 3.4 and 3.5 imply clearly that HA^m < C for all V E V. This finishes the proof of 
(|27|). The proof of (|28|) is quite similar and is therefore omitted. 

We now prove Lemmas p.3(-|3.5| to conclude the proof of Lemma 3.1. 



Proof of Lemma |3.3|. The Fourier transform of A| in k is given by 



C\{x, k) 



1 rdv{p) 



% 

Ai(a;, - + rfc, A: - ^, y) - Ai(a;, - + rfc, fc + |, V) 
£ 2 e 2 



The second term above may be written as 



i / dke-'''-y 



dV{p) 



;ip-{x/s+rk) 



/.^(ri.— / ...» 



dse"^ 



Ai(.T, ^+rk,k~^,V)^ Ai(.T, - + rfc, k + ?-,V) 
e 2 £ 2 



dVjq) 

{2TT)d ' 



,is{k—p /2)-q+i{x / e+rk) ■ q 



dke-'^-y 



dVjp) 
{2ttY ' 



^i'p-{x I e-\-rk) 



dVjq) 
{2Tr)<i ' 



,is{k-{'p/2)-q'\~i{x/6'\-rk)-q 



Kk + ^--^)-x{k + - + -) 



This is further transformed to 



dke-'^y 



dVjp) 
{2n)d ' 



,ip-{x/e-'rrk) 



dse'^ 



_^is{k-p/2) 
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dV{q) 
{2tiY'' 

(k-p/2-q/2) gia'-(fe-j9/2+q/2)| 



'i^/^+^f^y-i / dy'~X{x,y') 



^is(k+p/2) 



•9 y^y' 



.{k+p/2-q/2) _ iy'-{k+p/2+q/2)\ 



dVjp) 



,ip-(x/e) 



dse 



sQ 



dV{q) ix/e-q 

{2Tr)d 



A(x, y — rp — (r + s)q) 



X |_g-«sp-9/2-i(p+g)-(a-''p-(i-+s)«)/2 _|_ ^~isp-q/2-i(p~q)-{y-rp-(r+s)q)/2 
_j_ ^isp-q/2+i(p-q)-(y-rp-(r+s)q)/2 _ gisp-g/2+i(p+g)-(j/-rp-(r+s)g)/2 | 



Therefore we obtain 



\XUx,y,V)\<C I dre-"'-sup 

V 



\dVip)\ 
(27r)'' 



dse sup 



Vi 



\dViiq)\ 

{2ttY 



X \X{x,y~rp^{r + s)q)\+C / dre-"n'CA(x, y)! 



which is (I 



Proof of Lemma 3.4 



The first inequality in (bq) follows form the definition of | 



and 



the second is shown as follows. Let us define 

dp A . fc^ 



g{x, k) = LX{x, k) = 



Taking the Fourier transform in x and k we obtain 



-,k~p) [X{x,p) - X{x,k)] . 



dxdkdpdy' fyfk'^ -P^ 



{2n) 



2d 



■yR{- 



,k-p)X{x,y') 
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Integrating x out we obtain 



dkdpdy' -^k■yJ,f^P'^ 



We make a change of variables k' — k ~ p, p' = {k + p)/2 and drop the primes to get 



dkdpdy' _ 



2d 



• p, k)\{£„ y' ) 



^-ik-y' /2 _ ^ik-y'2 



dkdpdy' ds 



2d 



-isk-p—ip-y — ik-y/2-^ip-y 



xRis,k)X{^,y) 



g-ik-y' /2 _ ^ik-y j2 



We may now integrate p and y' out to obtain 

dkds 



Therefore we have 



<2 / ^|i?(s,fc)||A(e,y + sfc)|. 



rfCd2/|ff(^,2/)l <2||i?(s,p)|Ui / didym,y)\ 



and thus fB6f) holds 



Proof of Lemma 3.5. Clearly we have 



and thus it suffices to look at \y\ > 2. We observe that 

\dVip)\ 



Ir>s{x,y) = / dre sup 

-ISiy) V 
<^e-^S[y)^\\h 



{27rr Jo 



dse sup 

Vi 



(40) 



M^i(9)l,w 

L-^l^\X(a:,y~rp-ir + s)q)\ 



Therefore we have 

J dy Snp Ir>six,y) < Ca\\\\\L° 

Now we look at Ir<s ^-nd split it as: 



siy) 



dre sup 



\dV{p)\ 
(2^)^ Jo 



dse sup 



\^\Xix,y-rp-ir + sM 



Ir<s{x,y) 

Observe that 
so that 



It remains to bound /^^f . Note that for r,s < S{y), \p\, \q\ < L and A in the Schwartz class we 
have 



<i:.fsi^^y) + ir<li^^y)- 

dy sup I {x,y) < C„||A||lo 



\X{x,y — rp — {r + s)q)\ < sup |A(a;,z)|< sup |A(a;,z)|< 



\z-y\<2{r+s)L 



kl>|y|/2 
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Therefore we have 

y|>2 



dy 

\V\>2 Jo 



S{y) 



dre ™' sup 



S{y) 



dse sup 

Vi 



\dVi{q)\ 



X |A(a;, y — rp — {r + s)q)\ < 



\y\>2 \y 



<Cx. 



This finishes the proof of Lemma 3.5. 



3.3 The tightness of the measures Ve- 

The process We{t) generates a probabihty measure on the space C{[Q,T]; X) with the space 
X defined as before X ^ {W e S' : \\W\\a' < C}. This family is tight. 

Lemma 3.6 The family of measures is weakly compact. 

Proof. We follow the corresponding proof of Blankenship and Papanicolaou Q| for oscillatory 
ordinary differential equations with random coefficients. A theorem of Mitoma and Fouque 
IitI , implies that in order to verify tightness of the family it is enough to check that 
for each A e Ci([0, T], 5(M'* x M'')) the family of measures Ve on C([0,T];M) generated by the 
random processes W^{t) — {We{t),X) is tight. Tightness oiVe would follow from the following 
two conditions. First, a Kolmogorov moment condition |^ in the form 

E^^ {\{W,X)it) - {W,X){h)r \{W,X){h) -{W,X){sr}<Cx{t-s)'+^, 0<s<t<T (41) 

should hold with 7 > 0, /? > and C\ independent of e. Second, we should have 

lim limsupProb^- i sup \{W,X){t)\ > r \ = 0. 



R^oc e^O 



0<t<T 



The second condition holds automatically in our case since the process Wf^{t) is uniformly 
bounded for alH > and e > 0. In order to verify (|4l|), note that we have 

{W{t),X)^Gl^{t)-^e{W,X\)~e{W,Xl)+ f' ds{W,^ + k-V^X+CX){s) + V^ f' ds{W,C){s). 



The uniform bound (|2^) on (^^ and the bounds on ||A^ 2(011-4 in Lemma 3.1 imply that it suffices 
to check (^) for 

xe{t) = Gl{t) + ds{W, — + fc . V,A + CX){s). 



We have 



Ei^\xe{t)-xAs)\'\Ts] <2e\^ dr(M^,|^ + fc- V,A + £A)(t) t}^ 

+2E[\Gl^{t)-Gl^{s)\^\Ts]<C{t-sf + 2E{{Gl^) (t) - {G^ ){s)\Ts} . 
Here {G\ ) is the increasing process associated with G| . We will now compute it explicitly. 



First we obtain that 

d „p 
— E^' - 
dh wy,t 



^2{W,X,){W,^ + k-W^X, + ^IC[V,x/e]X,) + -Q[{W,Xe)^] 
h=0 £ 
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so that 



{W, X.fit) - [2{W, K) {s){W, -g^ + k- V,A, + ^IC[V, x/e]K) (s) + -Q [{W, \,)^] (s)] ds 



is a martingale. Therefore we have 



(GUt)) = / ds 



-Q[{W,X,)^]--{W,X,){W,QX, 

e £ 



with 



* ds {Q [{W, Xl)^] - {W, XI) {W, QXl){s)) +^ f dsH,{s) 
Jo 



H, = 2y^iQ[{W, XI) {W, XI)] - {W, XI) {W, QXl) - {W, Xl){W, QX\)) 
+ e {Q[{W, Xl)^] - 2{W, XI) {W, QA|)) . 



Lemma 3.1 and the boundedness of Q on L°°(V) imply that |iJe(s)| < C for all V G V. This 
yields 

E{{Gi:){t)-{Gl^){s)\Ts]<C{t~s) 

and hence 

E[\x,{t)-xM\^\^s} < Cit-s). 
In order to obtain ( ^l| ) we note that 

E^' {\xS) - xSi)nxSi) - ooeis)^} 

= E^' {E^^ { \x,{t) - xSl)r\^tA \Xe{tl) - X,(s)r} 

I \x,{t) - X,itl)f J^t,\\ \Xe{tl) - X,is)\^' 

< C{t - h)'>/^EP' {\xSi) - xeis)^} < C{t - h)^'^EP^ {EP^ {\x,{ti) - Xe{s)V \Ts}] 



<E' 



E' 



< C{t~ti)^/^EP^ 

< Cit-s)'-'. 



E' 



[\xSi) - x,{s)\^\Ts}] \<C{t~ ti)^^Hti - s)-*'^ 



Choosing now 7 > 1 we get ( |41| ) which finishes the proof of Lemma 3^. This also finishes the 



proof of Theorem 2.1 



4 Conclusions 

We have presented a proof of the transport limit (H) for the Schrodinger equation with a time- 
dependent random potential. Our proof is relatively simple and does not involve infinite Neu- 
mann expansions because it relies on the Markovian property of the potential, which allows us to 
construct approximate martingales and to show weak compactness of the family of probability 
measures generated by the dynamics ( p^ ) of the Wigner transform on C([0, T];A'). However, 
we only show convergence for the average Wigner distribution, which is the the first moment 
of Pe- We do not have a rigorous convergence result for the higher moments of the Wigner 
distribution and thus are not able to fully characterize the set of accumulation points of the 
family P^, although we believe that the limit measure P is unique, based on the formal analysis 
in [p8| of a similar problem in the white noise limit. 
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A The formal perturbation expansion 

We present the formal derivation of the transport equation for the hmit Wigner distribution 
W{t,x,k) similar to the one in |Q for a time-independent potential. Recall that the Cauchy 
problem for the Wigner distribution is 



dt 



1 f dV{t/e,p) 



1 



We{t,X,k-^)-We{t,X,k + ^)\ (42) 



WeiO,X,k) = W^{x,k). 



Here is the Wigner distribution of the family ^"(x), the initial data for (^). We will con- 
struct a formal perturbation expansion for We and derive the transport equation for the average 
W{t, X, k). We seek an expansion of with multiple scales 

Weit, x, k) = W(°)(t, x, k) -f y^W^^Ht, -,x,-,k)+ eW'^^\t, -,x,-,k) + ... (43) 

e e e e 

and introduce the fast time and spatial variables r = t/e, z — xje. We assume that the leading 
term W^") is deterministic and independent of the fast scale variables. We insert ( p3| ) into (^) 
and obtain at the order OiXj \fe): 

Then W^^^ has the form 

1 f dpduje'P-^+"^' V{uj,p) 



W^'^\t,T,X,Z,k) 



(27r)'*+i iuj + ip - k + S 



W^°\x,k^^)-W'-°\x,k + ^^. (44) 



Here (5 <C 1 is a regularization parameter that we will send to zero at the end of the calculation 
and V denotes the Fourier transform in time: 



V{uj,p) = / dse'''^'V{s,p). 
The term of order 0(1) in (nsl) gives 



dt 



d7 



-k-\/^W^°^ +k-S/^W^^'> 



) _ 1 /■ dpV{T,p)e^P-' 



W^(0)(fc_|)_^(0)(fc + |) 



We average the above equation assuming formally that E 
gives an equation for the leading order term W'^^^: 



aw(2 



{ d7 



k ■ 



(2) 



0. This 



dt 



1 f dpV^Typyp^ 



W^^\x,z,k-i^)-W^^\x,z,k + ^^ \. (45) 



The average on the right side of (^) may be computed explicitly using (^) and spatial homo- 
geneity (||): 



E- 



1 r dpV{T,p)e'P-' 
~i J {2^d 

dp ~ p^ — fc-^ 



{2tt 



W'^^\x,k-^)-W^^\x,z,k + ?^) 



,p-k) fw(")(p)- W(°)(fc) 
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as 5 ^ 0. Here R is the Fourier transform of R in time: 

R{uj,p) = / dse-"^'R{s,p). 
Js. 

Therefore we obtain the transport equation for the leading order term W^^^ : 



dt 



+ k ■ V^VF(°) 



i2n) 



,p^k) [Ty(")(p)-W^(o)(fc)' 



(46) 



The formal asymptotic expansion may not be justified but the final equation (^) for 
the expectation of the limit Wigner distribution E \W {t, x, k)} is correct. Moreover, the test 
functions that we used in our proof of Theorem |2.l| are based on the formal expressions for 
and W^^\ The role of the regularization parameter S is played by the spectral gap of the 
generator Q because of the bound ([T^. 
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